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Abstract 

The production and semi-leptonic decay of heavy quarks have been studied in the photoproduc- 
tion process e+p ^ e+ + dijet + + X with the ZEUS detector at HERA using an 
integrated luminosity of 38.5 pb"^. Events with photon-proton centre-of-mass energies, W^p, 
between 134 and 269 GeV and a photon virtuahty, Q^, less than 1 GeV^ were selected requir- 
ing at least two jets of transverse energy i?;^*^'-^-' > 7(6) GeV and an electron in the final 
state. The electrons were identified by employing the ionisation energy loss measurement. The 
contribution of beauty quarks was determined using the transverse momentum of the electron 
relative to the axis of the closest jet, The data, after background subtraction, were fit with 
a Monte Carlo simulation including beauty and charm decays. The measured beauty cross sec- 
tion was extrapolated to the parton level with the b quark restricted to the region of transverse 
momentum > p™™ = 5 GeV and pseudorapidity jry^j < 2. The extrapolated cross section 
is 1.6 ± 0.4 {stat.)~!io'^isyst.)~!iQil{ext.) nb. The result is compared to a perturbative QCD 
calculation performed to next-to-leading order. 
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1 Introduction 



High-energy collisions between a quasi-real photon, emitted by an incoming positron, and a proton 
can lead to the production of heavy quarks. Such processes allow a test of perturbative QCD 
(pQCD) since the mass of the heavy quark provides a hard scale. 

Measurements of charm production in 7p collisions at HERA have been made |||, |2| by reconstructing 
D*±(2010) mesons. The cross sections generally lie above the pQCD predictions. The study of 
beauty production is important since the heavier 6-quark mass provides a harder scale, thus making 
pQCD calculations more reliable. However, the higher mass and the smaller electric charge of the 
6-quark lead to cross sections in ep collisions that are typically two orders of magnitude smaller 
than that of charm. The first measurement of the cross section for beauty photoproduction has 
been performed using events with a muon and two jets Q and is higher than pQCD expectations. 
In pp interactions the measured beauty cross sections [Q,^, ^ significantly exceed the predictions. 

This paper presents a measurement of beauty production in photon-proton collisions using a sample 
of events each containing two jets and a candidate for the electron from a semi-leptonic decay of a 
heavy quark: 

e+p ^ e+(7) p ^ e+ +dijet + e" +X (1) 

Reaction (||) is isolated by statistically subtracting the hadronic background from an electron- 
enriched sample of events selected using measurements of the ionisation energy loss of charged 
particles. 



2 Experimental conditions 

The data used were collected by the ZEUS experiment during the 1996 and 1997 running pe- 
riods, when HERA operated with protons of energy Ep = 820 GeV and positrons of energy 
Eg = 27.5 GeV. The data for this study correspond to an integrated luminosity of 38.5 it 0.6 pb~^. 
A detailed description of the ZEUS detector can be found elsewhere @, . A brief outline of the 
components that are most relevant for this analysis is given below. 

Charged particles are measured in the central tracking detector (CTD) @, which operates in a 
magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists of 72 cylindrical 
drift chamber layers, organised in 9 superlayers covering the polar angle|^ region 15° < 9 < 164°. The 
transverse-momentum resolution for full-length tracks is Gp^/pT = 0.0058pT © 0.0065 © 0.0014/pr 
with pt in GeV. The energy loss of charged particles per unit track length, dE/dx, is also measured 



in the CTD |0| 



The high-resolution uranium-scintillator calorimeter (CAL) [p4| consists of three parts: the forward 
(FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part is subdivided transversely 
into towers and longitudinally into one electromagnetic (EMC) and either one (in RCAL) or two 
(in BCAL and FCAL) hadronic sections (HAG). The smallest subdivision of the calorimeter is 
called a cell. The electromagnetic section of the BCAL (BEMC) consists of cells of ~20 cm length 
azimuthally and mean width 5.5 cm in the Z direction at a mean radius of ~1.3 m from the beam 

^The ZEUS co-ordinates form a right-handed system with positive-^ in the proton beam direction and a horizontal 
X-axis pointing towards the centre of HERA. The nominal interaction point is at X = y = Z = 0. 
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line. These cells have a projective geometry as viewed from the interaction point. The CAL relative 
energy resolutions, as measured under test beam conditions are 18%/V^ for electrons and 35%/\/B 
for hadrons {E in GeV). Energy deposits in the CAL were used to measure the transverse energy 
and direction of jets. Cell clusters were also formed, which were then used to aid in the identification 
of electrons for the cross-section measurements. 

The luminosity was measured from the rate of the bremsstrahlung process e~^p — > e'^jp, where the 



photon was measured in a lead-scintillator calorimeter |12| located at Z = —107 m. 

A three- level trigger system was used to select events online [^, 13 1. At the third level, a cone 
algorithm was applied to the calorimeter cells and jets were reconstructed using the deposited 
energy and positions of the CAL cells. Events with at least two jets, each of which satisfied the 
requirements that the transverse energy exceeded 4 GeV and pseudorapidity was less than 2.5, were 
accepted. 



3 Analysis 

3.1 Offline cuts and event selection 



To suppress backgrounds from beam-gas interactions, cosmic rays and from deep inelastic scattering, 
the following cuts were applied: 



neutral current deep inelastic scattering events with a well measured scattered positron can- 



didate in the CAL were removed by cutting on the inelasticity, y, |14] which is the fraction of 



the electron energy carried by the photon in the proton rest frame. For an incoming positron 

2Ee 



of energy E^, y is estimated from ye = 1 — (1 — cos^g) where E'^ and 9'^ are the energy 



and polar angle of the outgoing positron. 

the requirement 0.2 < yjB < 0.8 was imposed, where yjB is the estimator of y measured 



from the CAL energy deposits according to the Jacquet-Blondel method |15|. The cut was 
imposed after correction for energy loss due to inactive material in the detector. This range 
in yjB corresponds to a photon-proton centre-of-mass energy, W^p from 134 to 269 GeV. 



These cuts restrict the photon virtuality, Q^, to less than 1 GeV^. The median value of about 
10~^ GeV^ was estimated from a Monte Carlo (MC) simulation. 

Jets were reconstructed using the KTCLUS [^] finder in its "inclusive" mode [|l^]. KTCLUS is a 
clustering algorithm which combines objects with small relative transverse energy into jets. The 
objects input to the jet algorithm may be hadrons in a simulated hadronic final state, the final- 
state partons of a pQCD calculation, or energy deposits in the CAL. After applying the KTCLUS 
jet algorithm to the calorimeter cells, the jet transverse energy was corrected for energy loss due 
to inactive material in the detector as a function of -^T CAL described in a previous 

ZEUS publication |jl^. After these corrections, all jets with \Tf^^\ < 2.4 and a transverse energy 
Sj?* > 6 GeV were kept. Each event was required to have at least two jets satisfying these criteria, 
with at least one jet having Ej^^ of more than 7 GeV. 
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3.2 Electron finding 



Electrons were identified in the CTD from the dE/dx of charged tracks, using a method of statistical 
subtraction |18, |l^]. Two samples of tracks were defined using information from the CAL: the first 
enriched with electrons with a background of hadrons (electron-enriched sample) and the second 
containing only hadrons (hadronic sample). The two samples were obtained by considering clusters 
in the calorimeter and performing a selection on the basis of the energies deposited in the EMC and 
HAC sections of the CAL. The CAL clusters matched to tracks were required to have a distance of 
closest approach, to the track of less than 20 cm. 



3.2.1 Electron-enriched and hadronic samples 

Using the matched track-cluster pairs, the electron-enriched and hadronic samples were defined by 
a cut on the ratio of electromagnetic energy, -BemCj to total energy, -Btot, of the clusters. The 
electron-enriched sample^ was required to have -EemcZ-E-tot > 0.9, while the hadronic sample was 
required to have -EemcZ-S'tot < 0.4, with a further requirement on the energy deposited in the 
hadronic section of the calorimeter, -Ehac > 0.3 GeV. These requirements on the hadronic sample 
efficiently rejected electrons; the residual contamination from electrons was estimated from photons 
converting into electron-positron pairs, reconstructed as described below, to be less than 0.03%. 
The selection criteria for the electron-enriched sample were 75% efficient for tagging electrons, as 
determined from the same sample of photon conversions. 

3.2.2 Measurement of dE/dx 

Charged particles traversing the CTD lose energy primarily by ionising gas in the detector. In 
order to estimate dE/dx for each track, the truncated mean of the anode- wire pulse heights was 
calculated, which removes the lowest 10% and at least the highest 30% depending on the number of 
saturated hits. For electrons traversing all superlayers, implying a maximum of 72 possible hits, on 
average 32 pulse-height measurements were retained. 

Since the CTD operates at ambient atmospheric pressure, the dE/dx calibration was changing 
throughout the measurement period. The measured dE/dx values were corrected by nor- 
malising to the average dE/dx for tracks around the region of minimum ionisation for pions, 
0.3 < p^^^ < 0.4 GeV, where the separation from other types of particles is particularly good. 
Henceforth dE/dx is quoted in units of mips - minimum ionising particles. 

The measured dE/dx value also depends on the polar angle. There is a trivial 1/cos^ dependence 
due to the path length which is corrected for in the subsequent plots. In addition, there is a 
dependence on 6 arising from the reduction in gain that occurs through screening of ions in the 
avalanche. This effect was studied using a sample of electrons originating from photons which 
converted in the beam-pipe via the 7 — > e~^e~ process. The candidate tracks were initially selected 

^Since the value of dE/dx is particle dependent, the relative fractions of pions, kaons and protons determine the 
background in the electron- or positron-enriched sample. This background can only be estimated using the hadronic 
sample. Monte Carlo studies showed that, for positively charged particles, the n/{K,p) ratio of the hadronic sample 
was markedly different to that of the hadronic background in the positron-enriched sample. This effect is caused by 
the differing cross sections for positive, compared to negative, low-energy pions, kaons and protons interacting with 
nuclei Therefore, since a reliable statistical subtraction was not possible, positrons were not considered further in 
this analysis. 
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on the basis of their distance of closest approach, vertex position and invariant mass. The quahty 

factor, D = {llxy / a^yY + {^G/^eY-, was calculated, where Axy is the separation of the tracks in 
two dimensions at the point at which their tangents are parallel, is the difference in polar angles 
and axy and cjq are the respective resolutions . The distribution of the quantity D is shown in 
Fig. |l]a, which demonstrates that real conversions, with net charge zero, tend to have lower values of 
D and as shown in Fig. ||b, low invariant masses of the electron-positron pair, Mg+g- . To achieve 
a pure sample of electrons, relatively hard cuts oi D < b and Mg+g- < 0.025 GeV were applied. 
Electron-positron candidates of net-zero and net-two units of charge were considered and the two 
subtracted as shown in Fig. Qa. The dE/dx distribution for the sample of clean electrons has a 
roughly Gaussian shape centred about dE/dx ~ 1.4 mips with width 0.14 mips, corresponding to 
a resolution of ~ 10%. Figure shows that the dE/dx value exhibits a strong dependence on the 
polar angle, 6^^'^, as expected from the space-charge effect [^]. The value of dE/dx at 9^^'^ = 7r/2 
is about 10% lower than the most forward and backward values in the range of 9^"^^ considered. 
Using this sample of electrons, correction factors were obtained which depend on 9^^^^. The dE/dx 
for electrons was corrected such that the mean was 1.4 mips. 



3.2.3 Tracking requirements 

Negatively charged tracks with a transverse momentum relative to the ep beam axis, Pj^, greater 
than 1.6 GeV were selected. Electron candidates were restricted to the central region of the detector, 
jj^trkj ^ -^ -^^ corresponding to 0.64 < 9^^^ < 2.50 radians, where the resolution in dE/dx is constant 
to within 10%. A small slice in rj^"^^ was removed from the analysis, corresponding to the region in 
which tracks were matched to clusters where the BCAL and RCAL meet, where CAL clusters were 
not well reconstructed. Figure ^ shows that for 1.6 < pj'^ < 10 GeV, all hadrons have an average 
dE/dx value well below that of electrons, thus allowing the separation of electrons from hadrons. 



3.2.4 Background electrons from converting photons 

The major source of background to electrons from semi-leptonic decays of heavy quarks (prompt 
electrons) comes from photon conversions in the detector. Electron-positron pairs were initially 
selected as discussed in Section p. 2. 2 . A loose cut of D < 15 was chosen. 



Tracks above a momentum of 0.2 GeV were reconstructed with an efficiency of greater than 95%. The 
number of potential conversion candidates not identified after this cut because the electron-positron 



pair was asymmetric in momentum was determined from a calculation of pair production |22] rather 
than by relying on a MC model. After the p^^^ > 0.2 GeV cut on the positron candidate was 
imposed, good agreement was seen between data and expectation for the shape of E^- /E^^ the 
fraction of the photon's energy carried by the electron, as shown in Fig. |^. This demonstrates that 
the identified electron-positron candidates originated from photon conversions. The requirement 
for the positron's momentum, p^^ > 0.2 GeV, was then removed in the calculation and the ratio, 
e{E^) = 7V(pfr > 1.6 GeV, p^"^ > 0.2 GeV)/7V(p|r > 1.6 GeV) of the two theoretical predictions 
(the solid and dashed curves in Fig. |^ determined, where is the number of photon conversions. 
The efficiency, e, averaged over the photon energy was found to be 83%. 

The overall efficiency for tagging the background due to electrons from photon conversions was 
determined from MC simulations. This was achieved by performing the analysis procedure with 
an inclusive MC conversion sample but demanding no prompt electron. The number of electrons 
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found by the analysis procedure, together with the number of electrons identified as coming from 
converting photons by the above two steps, was determined. The ratio of the number of identified 
conversion candidates to the total number of candidates passing the analysis cuts was 90 it 3% 
which in combination with e = 83%, leads to an overall efficiency for the identification of conversion 
candidates of 75 ± 3%. 

3.2.5 Electrons from Dalitz decays, 7r° — ^ e+e~7 

A substantial background is produced by Dalitz decays, vr'' e~^e~'y. To remove this background, 
all tracks in the electron-enriched sample that were not identified as coming from conversions, were 
combined with tracks of positive charge and the invariant mass formed. The positively charged 
tracks were required to have originated from the primary vertex, to have passed through at least 
the first three super layers of the CTD and to have pj}^ > 0.1 GeV. Since MC simulations showed 
that no fully reconstructed Dalitz decays survived a cut of Mg+g- > 0.2 GeV, an electron candidate 
was removed if any combination failed this cut. The overall rejection efficiency was 84 it 2%. 

3.2.6 Electron signal 

In order to subtract statistically the hadronic background from the electron-enriched sample, the 
dependence of the measured dE/dx on momentum and polar angle must be taken into account. 
As the p and 9 distributions for the electron-enriched and hadronic samples differ, subtraction of 
the hadronic background was carried out by reweighting the hadronic sample to give the same 
distribution as the hadronic background in the electron-enriched sample. This was achieved by 
binning the distributions in transverse momentum and pseudorapidity and normalising the hadronic 
sample to the electron-enriched sample in the region 0.5 < dE/dx < 1.1 mips. In this analysis, 64 
bins, with eight divisions in pj<^ and eight in rj^'^^, were used. 

Figure ^ shows the dE/dx distributions for both the electron-enriched and hadronic samples. The 
subtracted distribution together with the background from the photon conversions are shown in 
Fig. |5|b. Similar distributions are seen in the individual bins of pj'^ and r/*'"'^. In Fig. ^a, the two 
distributions have the same shape for dE/dx values below ~ 1.1 mips, but there is a clear excess 
of the electron-enriched sample over the hadronic sample at larger values of dE/dx, indicating the 
presence of electrons. The excess over the conversion signal was used to extract cross sections for 
electrons from heavy-quark decays. A cut at the mean value for electrons of dE/dx = 1.4 mips 
was made, which, assuming a symmetrical Gaussian distribution (as suggested by Fig. |2|a), has 
an efficiency of 50%. Varying the cut in steps of 0.01 mips between 1.3 and 1.5 mips gave the 
same results within the statistical uncertainties. Cutting at lower values of dE/dx gave no overall 
improvement in the statistical error. 

The number of identified electrons after this procedure is 1480 it 63, of which 537 it 29 were 
attributed to photons converting into electron-positron pairs, resulting in 943 it 69 electrons used 
for the cross section determination. The same statistical subtraction procedure was used in each 
bin to extract the differential cross sections. 
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4 Monte Carlo event simulation 



The acceptance and the effects of detector response were determined using samples of MC events. 
The programs Herwig 5.9 and Pythia 5.7 |2^, which implement the leading-order matrix 
elements followed by parton showers, were used. For all generated events, the ZEUS detector 
response was simulated in detail using a program based on GEANT3.13 |p5| . 

At leading order (LO), two types of processes can be distinguished: direct photon processes, in 
which the photon couples directly to a parton in the proton; and resolved photon processes, where 
the photon acts as a source of partons, one of which participates in the hard interaction. Samples 
of direct and resolved photon events, including heavy-quark excitation processes, were generated 



separately. For acceptance corrections, the MC events were generated with the CTEQ-4D |2(:] 
structure function for the proton and GRV-LO for the photon. The default quark masses were 
used in both Herwig and Pythia. 

For the fitting procedure and MC predictions, the GRV94-LO structure function for the proton 
and the GRV-LO for the photon were used. The quark masses were set to the nominal values of 
rriiy = 4.75 GeV and nic = 1-5 GeV. Samples with different input parton density functions and 
different quark masses were used to evaluate systematic effects. 

Both the shape and normalisation of the MC samples were compared to the differential cross sections. 
Fits to the data yielded the fraction of resolved photon processes as well as of beauty production. 
The cross-section value predicted by the MC model was also compared to the measured 6 — > e~ 
cross section. The MC predictions were then used to extrapolate the measured cross section to the 
parton level to facilitate a comparison with an NLO calculation. 



5 Event characteristics 

Comparisons of the distributions of kinematic quantities between the data and the Herwig and 
Pythia MC simulations, which include production of all five flavours of quark, are shown in Fig. ^. 
Both MC models describe the data reasonably well. 

The number of jets is more accurately described by the Pythia MC but the jet quantities themselves 
are well described by both MC programs. The general shape of the pj^^ distribution is well described, 
although there are slightly more events in the data at high p^^. The matching of the track-cluster 
pairs is also well described, as demonstrated by the distribution of d, the distance of closest approach 
of the track and the cluster. The separation of the jet and electron candidate in r/ — space. Re- jet, 
is peaked at low values with a flat tail at i?e-jet > Ij amounting to 10% of the sample; it is 
described reasonably by the MC programs. Confidence in the use of the MC programs for acceptance 
corrections was thus justified. 

Also shown in Fig. ^ is the contribution in the MC model from the semi-leptonic decays of beauty 
quarks. These events, in general, have jets with larger transverse energy and have a tendency to 
be more forward (i.e. nearer to the proton beam direction) in r/J*^*. The prompt electron is more 
separated from the jet, being produced with a higher transverse momentum and also in a more 
forward direction. 
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6 Cross section measurements 

The signal for beauty decays can be seen in Fig. |7|, in which the p^^ distribution, where is the 
momentum of the electron transverse to the axis of the jet to which it is closest, is compared to the 
MC prediction. The data peak at low p'^^, consistent with predominantly semi-leptonic decays of 
charm quarks, with other contributions mostly from r and rj decays contributing less than 3%. At 
high pj^^, the data fall less steeply than the predictions for charm. The data are consistent with a 
significant contribution from b decays. 



6.1 Differential cross sections 

Differential cross sections were determined as a function of pj^^ and x^^, the fraction of the photon's 
momentum contributing to the production of the two highest transverse energy jets. The 

variable x°^^ is defined as 
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where the sum is over the two jets of highest transverse energy. Cross sections for reaction (|i|) were 
measured in a restricted kinematic region with p^ > 1.6 GeV and \t]^ \ < 1.1. Differential 
electron cross sections, da/dpj?^ and da/dx^^ , were determined in the region < 1 GeV^, 

0.2 < y < 0.8, requiring events with at least two jets with > 7(6) GeV and Ir/-''^*! < 2.4. 

For a given luminosity, C, the cross section, fij, in bin i is given by cTj = Nf°'^'^ / (C ■ A), where 
jycorr -g acccptance-corrected number of electrons in the bin i and A is the bin width. The 
acceptance correction-factors were obtained from MC simulations using a bin-by-bin method. The 
reference MC model was Herwig, with Pythia used as a systematic check. At low p^', the value 
of the correction factor was 2.3, decreasing with increasing values of p"^^ to a minimum of 0.8. For 
the cross section as a function of x'^'^, the correction factors were in the range 1.1 — 2.7, increasing 
with increasing x'^^. 

The measured differential cross sections, da/dp"^^ and da/dx'^^, where the e~ comes from the semi- 
leptonic decay of a heavy quark, are shown in Figs. |^ and |8| respectively. The contribution from 



Dalitz decays was removed via the procedure and with the efficiency estimated in Section p.2.5| . The 
data are compared to the expectations of the Herwig MC simulation, which was area normalised 
to the data for a comparison of shape; the scaling factor was 3.8. 

Figure |8| shows a peak at high x"^*^, consistent with direct photon processes. However, the tail at 
low 2;°^** cannot be explained by direct processes alone. The Herwig prediction of 35% resolved 
photon contribution (including flavour-excitation processes) and 65% direct gives good agreement 
with the data. Fitting the shape of the direct and resolved photon MC distributions to the data 
gave a resolved component of 28±5(siai.)% (x^/ndf = 1.5). It can therefore be concluded that LO 
MC models require resolved photon processes to describe heavy quark production, as was observed 
in the ZEUS result for charm Hill. 
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6.2 Beauty cross section 



The beauty cross section was extracted by fitting the distribution of the data to the sum of 
contributions from beauty and charm. In the fit, the relative fraction of beauty and charm in 
the MC simulation was varied and the fraction of beauty processes in the kinematic region was 
extracted by minimising x^- The quoted beauty cross section below includes only the direct semi- 
leptonic decay from a b quark to an electron since the cascade decay, 6 — > c — > e~ is included in 
the background expected from charm because it exhibits a pj^^ spectrum more similar to a charm 
than a beauty decay. 

The percentage of beauty production was determined to be 14.7 ± 3.8 {stat.) % (xV^df = 1.1). 
Using this, the cross section for beauty production in the restricted kinematic region: < 1 GeV^, 
0.2 < y < 0.8, with at least two jets, E^"^^^ > 7(6) GeV and Iryj'^*! < 2.4 and a prompt 
electron with pj. > 1.6 GeV and \r]'^ \ < 1.1, was found to be 



a,V;Ze++dijct+e-+x = 24.9 ± QAt'ri pb. (3) 

The predictions from the two MC models are 8 pb for Herwig and 18 pb for Pythia, using the 
parameter settings quoted in Section ^ for the fitting procedure. The large difference in the MC 
predictions comes mainly from the different default treatments of as, different scales and hadro- 
nisation and the use of massive matrix elements in Herwig [23| and massless in Pythia ||2^ for 



the generation of flavour-excitation processes. Despite the cross-section differences, the predicted 
fractions of beauty production are 16.2% in Herwig and 19.5% in Pythia, both in reasonable 
agreement with the data value of 14.7 it 3.8 (stat.) %. Using the same procedure for charm, the 
measured cross section was found to be in reasonable agreement with the ZEUS measurement of D* 
production 

6.3 Systematic uncertainties 

A detailed study of possible sources of systematic uncertainty was performed for the differential and 
total cross-section measurements by varying cuts on the quantities to check the acceptance and by 
using the alternative MC model to check the stability of the fit. The largest contributions to the 
systematic error for the differential cross sections da / dp'^^ and da / dx'^^ were due to the uncertainty 
in the CAL energy scale and the use of Pythia instead of Herwig to correct the data. In the case 
of the measured beauty total cross section, the significant errors were: 

• using Pythia instead of Herwig gave an uncertainty of —8%; 

• the range for yjB was varied by the resolution of 8%, leading to an uncertainty of lii%; 

• requiring i?c-jct < 2 rather than no requirement gave an uncertainty of —10%. 

• the normalisation range of the hadronic background to the electron-enriched sample was 
changed to 0.5 < dE/dx < 1.04 mips, changing the cross section by —12%. 

• varying the distance of closest approach of the track-cluster pair, d, by it 5 cm changed the 
cross section by !!l3%. 
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• varying the cuts for the hadronic sample on the fraction -Bemc /-^TOT between 0.3 and 0.5 and 
for the electron-enriched sample between -E'emcZ-^'TOT > 0.85 and 0.95 gave an uncertainty 
of 

• using the CTEQ-4L proton structure function yielded a systematic error of —4%, while using 
the GS96-LO [^] photon structure function yielded a systematic error of —5%. 

All systematic uncertainties on the differential cross sections, excluding the correlated uncertainties 
due to the luminosity and hadronic energy scale, were added in quadrature. For the total cross 
section, all systematic uncertainties except that due to the luminosity were added in quadrature. 



7 Extrapolation to parton level cross section 
7.1 Extrapolation procedure 

The measured beauty cross section was extrapolated to the parton level in a restricted range of the 
transverse momentum and pseudorapidity of the beauty quark using the two MC models. The region 
in pseudorapidity for the beauty quark was defined such that the acceptance was reasonably constant 
over the range considered. The MC simulation predicts that 95% of events used in the cross-section 
determination of Eq. (^) have a 6-quark of transverse momentum > 5 GeV. Accordingly, the 
measured cross section was extrapolated to the parton level for the region pj- > p™™ = 5 GeV, 
|ry''| < 2, < I GeV^ and 0.2 < y < 0.8. The extrapolation factor /ext is then given by 



fe: 



b->e- 

" e+p^e+ +dijet+e- +X 



MC 



The value obtained using the Herwig MC generator was /ext = 6.8. To correct to the full b cross 
section the branching ratio for the process of 6 ^ e~ (10.73 it 0.35) % was used. The major 
source of uncertainty on the extrapolation factor arises from the uncertainty in the hadronisation, 
quantified by the use of Pythia, to be -26%; Pythia uses the Lund model while Herwig uses 
the cluster model. The other significant sources of uncertainty arise from varying the space-like 
evolution parameters^ (—11%), varying the mass of the b quark from 4.5 GeV to 5.0 GeV (12*^) and 
the branching ratio (±3%). Applying /ext and the branching ratio to the measured cross section 
gives 

<f-.e+bx = l.Q±OMstat.)t'oi{syst.)tli{ext.) nb 
where the last error given is the estimate for the error arising from the extrapolation procedure. 



7.2 Comparison with NLO predictions 



The cross section is compared in Fig. y to a NLO QCD calculation |30]. The prediction comes 



from a fixed-order calculation in which b quarks are not active partons in the proton and photon 

^For the central analysis, the default setting, denoted by the variable ISPAC = 0, allows the backward evolution 
down to some cut-off, QSPAC, which was set to 1 GeV. By changing to ISPAC = 1, the backward evolution continues 
to the infra-red cut-off, but the parton density functions are frozen at the value QSPAC, which is again set to 1 GeV. 
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structure functions and are generated dynamically in the hard sub-process. The NLO calculation 
uses the MRST99 and GRV-G HO [^] parton-density parametrisations for the proton and 
photon, respectively. For the central prediction, the renormalisation and factorisation scales are set 



to the transverse mass, = fip = uit = y fn^ +Pt, where rrih = 4.75 GeV. The predicted cross 

section in the region > 5 GeV, < 2, < I GeV^ and 0.2 < y < 0.8 for the above 
settings is 0.64 nb. The values produced by variation of the mass of the beauty quark (4.5 and 
5.0 GeV) and the renormalisation and factorisation scales {2mT and mr/2) are shown in Fig. ^ as 
the upper and lower predictions. Using other sets of parton density functions, e.g. MRST99 {g |), 
MRST99 {g T), MRST99(as ii) and MRST99(as TT), in which the extremes of an allowed range 
of values are taken [^], resulted in variations of the NLO predictions for > 5 GeV of within 



lb 4%. Using CTEQ 5M1 |3^] and GRV98 HO [33| proton structure functions changes the prediction 
by +2% and —5%, respectively. 

The extrapolated cross section lies somewhat above the central NLO prediction, consistent with the 
general observation that NLO QCD calculations underestimate beauty production both in hadropro- 
duction IQ, H, |6| and photoproduction at HERA ||] . 



8 Summary 



The production and semi-leptonic decay of heavy quarks has been studied in the photoproduction 
process e~^p ^ e"*" + dijet + e~ + X with the ZEUS detector at HERA using an integrated 
luminosity of 38.5 pb~^. Events with photon-proton centre-of-mass energies, W^yp, in the range 
134 < W^p < 269 GeV and a photon virtuality, < I GeV^, were selected with at least 

two jets of transverse energy E^^^^^'^^ > 7(6) GeV and an electron in the final state. The beauty 
cross section was measured to be o^g+p^g+^^^^Q^_^^-_^_x = 

24.9 ± 6.4l^:^ pb. This cross section was 
extrapolated to the parton level with a b quark restricted to the region of transverse momentum 
Pj^>5 GeV and pseudorapidity |?7^| < 2 for the events with < 1 GeV^ and 0.2 < y < 0.8. 
The value obtained, crlp\,f,+i)X ~ ^-^ ^ ^-^ (stat.) ~^^'^{syst.)~^Ql{ext.) nb, lies somewhat above 
the NLO predictions, in agreement with results both from hadroproduction and photoproduction. 
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Figure 1: (a) The distribution of the quahty factor, D, (see Section [3. 2. 2[) for conversion candidates, 
(b) the invariant mass, Mg+g-, for conversion candidates. In (a) and (b), the conversion candidates 
resulting in pairs with zero net charge are shown as points; those pairs having non-zero net charge 
are shown as the crosses. 
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Figure 2: (a) The dE/dx distribution of photon-conversion candidates and (b) the dependence 
on the polar angle, 6^"^^, for electrons. In (a), photon-conversion candidates having two tracks of 
opposite charge are shown as open circles whereas those with tracks of the same sign are shown as 
the crosses; the solid circles show the difference between these two distributions. A Gaussian fit is 
shown in (a) for illustration. 
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Figure 3: The measured distribution of dE/dx against momentum, p*"^*^, for negative tracks in the 
range \if'^^\ < 1.1, as for the analysis. The curves show the expected average values for particular 
types of particles as derived from the Bethe-BIoch formula |34]. The events are a sub-sample of 
those that pass the dijet trigger requirements. 
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Figure 4: Comparison of the number of electrons found in the data (points) with the prediction 
from pair production (sohd hne) for conversions in which both tracks have a momentum greater 
than 0.2 GeV. The prediction for no cut on the momentum of the tracks is also shown as the dashed 
line. 
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Figure 5: (a) The dE/dx distribution for the hadronic sample (crosses) and electron-enriched (open 

squares) sample normalised to each other in the hatched region shown, (b) The difference between 
the electron-enriched and hadronic samples (solid circles), together with the background arising 
from photon conversions (open circles). The data in the region with dE/dx > 1.4 mips shown by 
the dashed line were used to extract the results in this paper. 
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Figure 6: Comparison of data (points) with Herwig (solid line) and Pythia (dashed line) MC 
expectations for: (a) jet transverse energy, -Ej?*; (b) pseudorapidity of the jet, if^^; (c) number 
of jets in the event, A'"j^*; (d) the Jacquet-Blondel estimator of y, j/jb; (e) the separation of the 
electron track and jet, i?e-jet; (f) the distance of closest approach of the track-cluster pair, d; (g) 
the transverse momentum of the track, pj^^, and (h) the pseudorapidity of the track, 77*'"'^. The MC 
prediction is area normalised to the number of data events. The contribution of events with an 
electron from the semi-leptonic decay of a 6 quark to the MC prediction is shown as the hatched 
histogram. The normalisation for the MC from 6-decays is the same as for the total sample and the 
fraction is that predicted in Herwig. 



18 



ZEUS 




12 3 4 



p;«'(GeV) 

Figure 7: The differential cross-section da/dp^r^^ for heavy quark decays. The inner error bars 
represent statistical errors and the outer bars statistical and systematic errors added in quadrature. 
The effect of the energy scale uncertainty is shown as the shaded band. The experimental data is 
compared to the Herwig MC prediction (solid line) for all quark flavours, which has been fitted 
to the data and scaled up by a factor of 3.8. The components from the beauty processes (forward- 
diagonally hatched histogram) and from charm (backward-diagonally hatched histogram) predicted 
by the Herwig MC model are indicated separately. 
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Figure 8: The differential cross-section da/dx'^^ for heavy quark decays. The inner error bars 
represent statistical errors and the outer bars statistical and systematic errors added in quadrature. 
The effect of the energy scale uncertainty is shown as the shaded band. The experimental data is 
compared to the HERWIG MC prediction (solid line) for all quark flavours, which has been fitted to 
the data and scaled up by a factor of 3.8. The direct (forward-diagonally hatched histogram) and 
resolved photon components (backward-diagonally hatched histogram) are indicated separately. 
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Figure 9: The extrapolated b cross section at a fixed p™™ value compared with theoretical predictions 

plotted as a function of p™™. The inner error bars represent the statistical error and the outer 
bars statistical, systematic and extrapolation errors added in quadrature. The curves show the 
predictions from NLO QCD for varying 6-quark mass and varying factorisation and renormalisation 

scale rriT = y + Pt- 
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